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1. Introduction 

The Riemann hypothesis (RH) is the statement that the complex zeros of the classical zeta 
function all have imaginary part equal to 1 /2. It was first suggested by Riemann in his famous 
memoir in 1859 [fl~|. The RH is important for its connection with the distribution of prime 
numbers ElO . The average number of primes less than a given number x, which is denoted as 
7i(x), behaves asymptotically as xj log x. This statement is called the Prime Number Theorem 
(PNT) and was proved independently by Hadamard and de La Vallee-Poussin in 1896 [[4j|5]. 
The truth of the RH implies that the fluctuations of n(x), around its average value, behaves 
asymptotically as x 1 ^ 2 log x, which also gives the best possible bound for the error of the PNT. 
The RH is not an isolated property of a particular function, but it holds for the Dirichlet L- 
functions, for curves over finite fields, etc. It is expected that a proof of the RH for the zeta 
function, will be generalizable to other L-functions as well. However the consensus is that 
some key idea is required for this goal. 

One of the most promising pathways for a proof of the RH was suggested by Polya and 
Hilbert around 1910, but never published apparently. The suggestion is that there exist a 
selfadjoint operator H, whose spectrum contains the imaginary part of the Riemann zeros. 
The selfadjointness of such an operator would inmediately prove the RH: 

If ^ + i£„J =0, Vn => H \ifr„) = E„ \if, n ) => E„ = E*„, Vn (1) 

In the latter equation one must exclude the trivial Riemann zeros, so a most appropiate 
formulation fo the problem is to use the function, defined as [4 | 

?(i) = i*(*-i)ff-' /2 r(!)f(*) (2) 
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whose zeros are those of ((s) except the trivial ones. Polya and Hilbert conjecture is known 
as the spectral approach to the RH, and, as we shall see later on, it is supported by several 
"phenomenological" results and heuristic arguments, which suggest that the operator H is 
the quantum Hamiltonian of a physical system. This is one of the reasons why the RH 
has attracted the interest of physicist working in disciplines apparently unrelated to Number 
Theory (see [6 1 for an extensive list of references concerning several approaches to the RH). 

Assuming the RH (i.e. E* = E„, Vn), one can define a diagonal operator H whose entries 
are E n , but nothing is learned from this construction. Eq.([T]i, implies that H must encode in 
itself, the zeta function £(1/2 + iE), without assuming the truth of the RH, which will be a 
consequence of its selfadjointness. If £(1/2 + iE), or £(1/2 + iE), where a polynomial in E, 
this encoding could be realized by a finite dimensional matrix H whose characteristic poly- 
nomial were proportional to £(1/2 + iE). The Euler product formula of the zeta function, in 
terms of the prime numbers, implies that H must also know about these numbers. Hence, the 
relation found by Riemann, between prime numbers and zeros of the zeta function, will be 
justified from the common dynamical origin of these quantities. The precise mathematical 
formulation of these relations is given by the, so called, trace formulas in Number Theory 
and Quantum Chaos [7 |. The simplicity of the definition of the zeta function, as the series 
£(*) = l/" 1 CR e s > 1)' has lead some researches, as Berry, to suggest that the Hamil- 
tonian H may have a subtle but simple definition, allowing the observation of the Riemann 
zeros as spectral lines in an experimental set up [8 |. The existence of such a "Riemann cal- 
culator" would place Number Theory in the realm of Quantum Mechanics with far reaching 
consequences. We shall briefly summarize below some phenomenologial and heuristic hints 
that support the spectral approach to the RH. 



Consider a compact Riemann surface with negative curvature. This surface can be con- 
structed as the complex upper plane divided by a discrete subgroup of the modular group 
PSL(2, R), and it is equipped with the Poincare metric. A classical problem is to determine 
the lengths i p of the primitive periodic orbits (p.p.o.), that is, the geodesies on this surface. A 
quantum problem is to find the spectrum of the Laplace-Beltrami operator A = y 1 (d\ + d 2 ), 



Selberg's trace formula establishes a relation between the momenta k„ and the length of 
the geodesies ( p 



2. Selberg's trace formula (1956) 



- A if/„(x, y) = E„ i// n (x, y), 




(3) 



CO 



g(n£ p ) 




2 smh(n{ p /2) 



(4) 



p.p.o. n=l 



where h(k) is a test function, g(k) its the Fourier transform and iu(D) is the area of the fun- 
damental domain D describing the Riemann surface. Selberg also defined a zeta function in 
terms of the lenghts ( p as 
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z(s) = nn (i -^" (5+m)) (5) 

p.p.o. m-0 

in close analogy to the Euler' s product formula of the zeta function J4): 

«*) = n y~~^' Re s > 1 (6) 

where the product is over all the prime numbers p. Selberg zeta function satisfies a RH which 
can be proved. The trivial zeros of Z(s) are s„ = —n (n = 0, 1, ... ), and the non trivial 
ones are s n — i + Since are real numbers, any complex zero of Z(s) lies on the line 
Re s — 1 /2. The functions g(s) and Z(s) both satisfy functional equations that relate their 
values to £(1 - s) and Z(l - s), respectively. Finally, Selberg' s trace formula is reminiscent 
to the Riemann- Weil explicit formula relating the prime numbers p and the imaginary part of 
the Riemann zeros y„ iflOl : 

(7) 

where the notations are as in eq.©. A comparison between eqs.(@]i and © suggests that 
prime numbers and primitive geodesies are in one-to-one correspondence, such that: l p <-> 
log p. This correspondence also underlines the Quantum Chaos approach to the RH reviewed 
later on. However, the analogy between these two formulas fails in two respects: i) the term 
1/(2 sinh(n£ p /2)), with the identification ( p — log p, only converges to p~"l 2 for large values 
of p, and ii) the factor of -2 in the last term of eq.(|7]i as compared to the factor one, in the 
last term of eq.©. The difference in signs of the two terms finds an explanation in Connes 
spectral realization of Riemann zeros (see below). Finally, we observe that the imaginary 
part of the Riemann zeros, y n , seem to correspond to the momenta k„ rather than to the 
eigenenergies E„. This suggest that the Riemann Hamiltonian is probably related to a first 
order linear operator, as we shall see in the discussion of the H — xp Hamiltonian. 

3. Random matrix theory and Quantum Chaos (70V80' s) 

In 1973 Montgomery, asuming the RH, proved that the imaginary part of the Riemann zeros, 
y n , where distributed at random, according the gaussian unitary ensemble distribution (GUE) 
of Random Matrix Theory (RMT) ifTTI . This result found strong numerical confirmation 
by Odlyzko in the 80' s, who computed trillions of Riemann zeros ( near the 10 12 -th zero 
and near the 10 20 - th zero) lfl2l . These phenomenological findings means that the statistical 
properties of the Riemann zeros are similar to those of the eigenvalues of large hermitean 
matrices, in particular the property of level repulsion. There are three universality classes 
of random matrices corresponding to orthogonal (GOE), hermitean (GUE) and symplectic 
matrices (GSE). The GUE statistics corresponds to random systems where time reversal is 
broken, which gives a strong indication that the Riemann Hamiltonian H must break this 
symmetry. 
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A further step along this direction was taken by Berry, who noticed a formal analogy 
between the fluctuations of the Riemann zeros and the fluctuations of the energy levels of 
quantum chaotic system around their average values [8]. The latter fluctuations are given by 
the semiclassical Gutwiller formula, 

1 v~i x~* sm(mETy) 

nqca(E) = -vv - . - (8) 

;r z—i z_i 2m smh(mAy 2) 

y 111=1 ' * 

where E is an eigenenergy, y is a primitive periodic orbit, T y its period and A 7 its Lyapunov 
exponent. The sum over m corresponds to the repetitions of the primitive orbits. The fluctua- 
tion part of the Riemann zeros is given by 

n *—l mp m ' z 
p m=l r 

where the sum is over the prime numbers p. Comparing ([8]l and ©, Berry conjectured 
the existence of a classical chaotic Hamiltonian whose primitive periodic orbits, y, would 
be labelled by the prime numbers p, with periods T p = log p, and instability exponents 
Hp = ± log p. Moreover, since each orbit is counted once, the Hamiltonian must break time 
reversal (otherwise there would be a factor 2/n, in front of eq. (O instead of 1/tt). This 
analogy is reminiscent to the one existing between Selberg trace formula and the Riemann- 
Weil formula, and also suffers from a "sign problem" and "asymptotic problem" as observed 
before. The connection with Quantum Chaos also explained some numerical discrepancies 
found by Odlyzko between RMT and the statistics of zeros for long range spectral correla- 
tions, which are due to the shortest periodic orbits, where universality no longer holds. They 
were explained by Berry , Keating and Bogomolny |[T3l[T4]| . All these results put on a more 
firm basis the Polya-Hilbert conjecture giving further clues on the structure of the dynamical 
system behind the Riemann zeros. 



4. The Hamiltonian H = xp (1999) 

In 1999 Berry and Keating, and Connes suggested that the Riemann zeros are related to the 
classical Hamiltonian H c i = xp, where x and p are the position and momenta of a particle 
movin in ID lfT31[T6ll . The classical trayectories of this Hamiltonian are hyperbolas in the 
phase space 



x(t) = x e', p(t) = po e 1 (10) 

and therefore unbounded, which would then imply a continuous spectrum rather than a dis- 
crete spectrum associated to the Riemann zeros. The connection with the latter arises in two 
possible ways depending on two different regularizations of the phase space. Berry and Keat- 
ing introduced a minimal length € x and a minimal momenta i p , whose product is the Planck 
quantum £ x ( p = 2nTi lfT5l . In terms of these quantities, they imposed \x\ > l x and \p\ > ( p , 
so that the trayectories are now bounded. The semiclassical number of states is given by the 
area below the hyperbola and above the boundaries |*| = l x and \p\ - ( p , and it is given, in 
units Ti — 1, by 
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(11) 



Rather surprisingly, this result coincides, asymptotically, with the average number of Rie- 
mann zeros up to a height E in the critical strip (i.e. < %s < 1,0 < 3 s < E). The 
constant in Riemann' s formula is actually 7/8, which can be obtained by taking into account 
a Maslow phase contributing -1/8 to eq. (fTTT l. due to the fact that the particle only travels 
one quadrant of the phase space. Unfortunately, eq. dTTb has remained so far heuristic since it 
is not supported by a quantum mechanical model (see later). 

Connes regularization is based on the restrictions \x\ < A and \p\ < A, where A is a cutoff, 
which is taken to infinity at the end of the calculation 1161 . The semiclassical number of 
states is computed as before yielding, 



The first term on the RHS of this formula diverges in the limit A — > oo, which corresponds 
to a continuum of states. The second term is minus the average number of Riemann zeros, 
which according to Connes, become missing spectral lines in the continuum. This is the, so 
called, "absortion" spectral interpretation of the Riemann zeros, as opposed to the standard 
"emission" spectral interpretation where they form a discrete spectrum. The minus sign in 
eq.(fT2l could also be related to the minus sign in the trace formulas discussed earlier. Connes 
interpretation has however two drawbacks. First of all, the average number of Riemann zeros 
is not fully obtained in eq (fT2l . The term EI2n\og l/2n, actually cancells between the first 
and second summands in this formula. Other objection is that the second term in (TTZb is 
simply a finite size correction of discrete energy levels, where no lines are missing, and the 
same remains true in the continuum limit. 

The H = xp model was modified in references ifFTl . adding a non local interaction 
suggested by a relation of this model to a BCS model of superconductivity with a cyclic 
renormalization group. The spectrum of this interacting xp-model is a continuum where the 
Riemann zeros are embedded as bound states. This result reconciles the Berry-Keating and 
Connes spectral interpretations. However, the non locality of the interaction implies that the 
Hamiltonian has no classical limit, and consequently its relation to classical chaotic dynami- 
cal systems remains unclear. Moreover the prime numbers do not appear in this construction, 
which as we saw in previous sections, is an important ingredient of the trace formulas. 

A different route to H — xp was suggested in reference [18| which will bring us to more 
familiar territores in Physics. 



Let us consider a charged particle moving in a plane under the action of a perpendicular mag- 
netic field and an electrostatic potential with a saddle shape [ 18 1. The Langrangian describing 
the dynamics is given, in the Landau gauge, by 




(12) 



Landau levels and Riemann zeros (2008) 



n A* 2 2 eD 

£ = r (x + y ) yx - eAxy 



(13) 
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where 11 is the mass, e the charge, B the magnetic field, c the speed of light and A a coupling 
constant parameterizing the electrostatic potential. There are two normal modes with real, 
u) c , and imaginary, w/,, angular frequencies, describing cyclotronic and a hyperbolic motions 
respectively. In the limit where lo c » u>h, only the Lowest Landau Level (LLL) is relevant 
and the effective Lagrangian becomes 

Hy lHc\ 1 ^ 2 Ac 



Xeff = px - \o)h\xp, p = ~p' { ^\e~B~) ' ~ ~B ^ 

where I is the magnetic length, which is proportional to the radius of the cyclotronic orbits in 
the LLL. The coordinates x and y, which commute in the 2D model, after the proyection to 
the LLL become canonical conjugate variables, and the effective Hamiltonian coincides with 
the xp Hamiltonian introduced by Berry, Keating and Connes, where the energy is measured 
in units of /z|w/,|. This realization of the xp Hamiltonian allow us to interpret the semiclas- 
sical quantization of these authors in the language of the Landau model. In particular, the 
semiclassical counting of states in the xp model follows from the counting of quantum fluxes 
in a certain area of the x — y plane. If the plane is infinite, then the number of states in 
the LLL will is also infinite. To have a finite number of states we put the particle into a 
box: |x| < L, \y\ < L, which reproduces Connes regularizations conditions. The number of 
semiclassical states with an energy between and E is given by 

-^^i^-ihi- 1 ) (15) 

which agrees with Connes eq. (fT2l . The classical energy is given by E — xy/{ 2 (in units of 
H\coi,\), and it attains its maximum value at E mm = L 2 /{ 2 . Plugging this value into (15[ yields 
«sm(£max) = L 1 j7.nl 1 , which is the number of quantum fluxes in the first quadrant. This 
semiclassical results can be derived from the quantization of the model. Indeed, the energy 
levels follows from the identification of the wave function at the boundaries x — L and y = L 
(up to a phase) 

Mx, L) = e'^MU x) ^YTW) \2fi) = 1 (16) 

Taking the logarithm on the RHS of (fTol l. one gets the smooth part of the Riemann formula, 
whose asymptotic expansion coincides with eq. Sl5[ . The 2D formulation of the xp model 
allows one to convert the hyperbolic orbits into periodic orbits by means of the boundary 
condition ( fTol l. It would be extremely interesting to derive Berry-Keating regularization in the 
Landau version of the model. To achieve this goal, one must inject the particle approaching 
the boundary at y = i y , back to the boudary at x = l x , so that the orbits become periodic. 
The Berry-Keating semiclassical arguments suggest that the spectrum will be discrete and 
associated to the smooth Riemann zeros. Preliminary results suggest that this possibility 
can indeed be realized. Of course the main problem that remains is the construction of the 
Hamiltonian giving rise to the exact Riemann zeros. It is not clear at the moment wether 
this can be done, but preliminary results suggest that the higher Landau levels may play a 
role IIT81 . In any case, the Landau model formulation of the H — xp Hamiltonian provides a 
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promissing new avenue where to explore the fascinating problem of a physical interpretation 
of the Riemann zeros, and perhaps a physicist proof of the RH. 
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